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MODIFIED INTERLEUKIN-1 RECEPTOR ANTAGONIST (IL-1RA) 

WITH REDUCED IMMUNOGENICS Y 



FIELD OF THE INVENTION 

5 The present invention relates to polypeptides to be administered especially to humans and 
in particular for therapeutic use. The polypeptides are modified polypeptides whereby the 
modification results in a reduced propensity for the polypeptide to elicit an immune 
response upon administration to the human subject. The invention in particular relates to 
the modification of human interleukin- 1 receptor antagonist (IL-1RA) to result in EL- IRA 

10 protein variants that are substantially non-immunogenic or less immunogenic than any 
non-modified counterpart when used in vivo. The invention relates furthermore to T-cell 
epitope peptides derived from said non-modified protein by means of which it is possible 
to create modified IL-1RA variants with reduced immunogenic ty. 

1 5 B ACKGRO UND OF THE INVENTION 

There are many instances whereby the efficacy of a therapeutic protein is limited by an 
unwanted immune reaction to the therapeutic protein. Several mouse monoclonal 
antibodies have shown promise as therapies in a number of human disease settings but in 
certain cases have failed due to the induction of significant degrees of a human anti- 

20 murine antibody (HAMA) response [Schroff, R. W. et al (1985) Cancer Res. 45: 879-885; 
Shawler, D.L. et al (1985) J. Immunol. 135: 1530-1535]. For monoclonal antibodies, a 
number of techniques have been developed in attempt to reduce the HAMA response 
[WO 89/09622; EP 0239400; EP 0438310; WO 91/06667]. These recombinant DNA 
approaches have generally reduced the mouse genetic information in the final antibody 

25 construct whilst increasing the human genetic information in the final construct. 

Notwithstanding, the resultant "humanized' 9 antibodies have, in several cases, still elicited 
an immune response in patients [Issacs LD. (1990) Sem. Immunol* 2: 449, 456; Rebello, 
P.R. et al (1999) Transplantation 68: 1417-1420]. 

30 Antibodies are not the only class of polypeptide molecule administered as a therapeutic 
agent against which an immune response may be mounted. Even proteins of human 
origin and with the same amino acid sequences as occur within humans can still induce an 
immune response in humans. Notable examples include the therapeutic use of 
granulocyte-macrophage colony stimulating factor [Wadhwa, M. et al (1999) Clin. 

1 
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Cancer Res. 5: 1353-1361] and interferon alpha 2 [Russo, D. et al (1996) Bri. J. Haem. 
94: 300-305; Stein, R. et al (1988) New Engl. J. Med. 318: 1409-1413] amongst others. 

A principal factor in the induction of an immune response is the presence within the 
5 protein of peptides that can stimulate the activity of T-cell via presentation on MHC class 
II molecules, so-called "T-cell epitopes. Such potential T-cell epitopes are commonly 
defined as any amino acid residue sequence with the ability to bind to MHC Class II 
molecules. Such T-cell epitopes can be measured to establish MHC binding. Implicitly, a 
"T-cell epitope" means an epitope which when bound to MHC molecules can be 
10 recognized by a T-cell receptor (TCR), and which can, at least in principle, cause the 
activation of these T-cells by engaging a TCR to promote a T-cell response. It is, 
however, usually understood that certain peptides which are found to bind to MHC Class 
II molecules may be retained in a protein sequence because such peptides are recognized 
as "self* within the organism into which the final protein is administered. 

15 

It is known, that certain of these T-cell epitope peptides can be released during the 
degradation of peptides, polypeptides or proteins within cells and subsequently be 
presented by molecules of the major histocompatability complex (MHC) in order to 
trigger the activation of T-cells. For peptides presented by MHC Class 11, such activation 
20 of T-cells can then give rise, for example, to an antibody response by direct stimulation of 
B-cells to produce such antibodies. 

MHC Class II molecules are a group of highly polymorphic proteins which play a central 
role in helper T-cetl selection and activation. The human leukocyte antigen group DR 

25 (HLA-DR) are the predominant isotype of this group of proteins and are the major focus 
of the present invention. However, isotypes HLA-DQ and HLA-DP perform similar 
functions, hence the resent invention is equally applicable to these. The MHC class II DR 
molecule is made of an alpha and a beta chain which insert at their C-termini through the 
cell membrane. Each hetero-dimer possesses a ligand binding domain which binds to 

30 peptides varying between 9 and 20 amino acids in length, although the binding groove 
can accommodate a maximum of 1 1 amino acids. The ligand binding domain is 
comprised of amino acids 1 to 85 of the alpha chain, and amino acids 1 to 94 of the beta 
chain. DQ molecules have recently been shown to have an homologous structure and the 
DP family proteins are also expected to be very similar. In humans approximately 70 
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different allotypes of the DR isotype are known, for DQ there are 30 different allotypes 
and for DP 47 different allotypes are known. Each individual bears two to four DR 
alleles, two DQ and two DP alleles. The structure of a number of DR molecules has been 
solved and such structures point to an open-ended peptide binding groove with a number 

5 of hydrophobic pockets which engage hydrophobic residues (pocket residues) of the 
peptide [Brown et al Nature (1993) 364: 33; Stern et al (1994) Nature 368: 215]. 
Polymorphism identifying the different allotypes of class II molecule contributes to a 
wide diversity of different binding surfaces for peptides within the peptide binding grove 
and at the population level ensures maximal flexibility with regard to the ability to 

10 recognize foreign proteins and mount an immune response to pathogenic organisms. 

There is a considerable amount of polymorphism within the ligand binding domain with 
distinct "families 11 within different geographical populations and ethnic groups. This 
polymorphism affects the binding characteristics of the peptide binding domain, thus 
different "families" of DR molecules will have specificities for peptides with different 

15 sequence properties, although there may be some overlap. This specificity determines 
recognition of Th-cell epitopes (Class II T-cell response) which are ultimately responsible 
for driving the antibody response to B-cell epitopes present on the same protein from 
which the Th-cell epitope is derived. Thus, the immune response to a protein in an 
individual is heavily influenced by T-cell epitope recognition which is a function of the 

20 peptide binding specificity of that individual's HLA-DR allotype. Therefore, in order to 
identify T-cell epitopes within a protein or peptide in the context of a global population, it 
is desirable to consider the binding properties of as diverse a set of HLA-DR allotypes as 
possible, thus covering as high a percentage of the world population as possible. 

25 An immune response to a therapeutic protein such as the protein which is object of this 
invention, proceeds via the MHC class II peptide presentation pathway. Here exogenous 
proteins are engulfed and processed for presentation in association with MHC class H 
molecules of the DR, DQ or DP type. MHC Class II molecules are expressed by 
professional antigen presenting cells (APCs), such as macrophages and dendritic cells 

30 amongst others. Engagement of a MHC class II peptide complex by a cognate T-cell 
receptor on the surface of the T-cell, together with the cross-binding of certain other co- 
receptors such as the CD4 molecule, can induce an activated state within the T-cell. 
Activation leads to the release of cytokines further activating other lymphocytes such as B 
cells to produce antibodies or activating T killer cells as a full cellular immune response. 

3 
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The ability of a peptide to bind a given MHC class II molecule for presentation on the 
surface of an APC is dependent on a number of factors most notably its primary 
sequence. This will influence both its propensity for proteolytic cleavage and also its 
affinity for binding within the peptide binding cleft of the MHC class II molecule. The 
5 MHC class II / peptide complex on the APC surface presents a binding face to a particular 
T-cell receptor (TCR) able to recognize determinants provided both by exposed residues 
of the peptide and the MHC class II molecule. 

In the art there are procedures for identifying synthetic peptides able to bind MHC class II 
10 molecules (e.g. W098/52976 and WO00/343 1 7). Such peptides may not function as T- 
cell epitopes in all situations, particularly, in vivo due to the processing pathways or other 
phenomena. T-cell epitope identification is the first step to epitope elimination. The 
identification and removal of potential T-cell epitopes from proteins has been previously 
disclosed. In the art methods have been provided to enable the detection of T-cell epitopes 
1 5 usually by computational means scanning for recognized sequence motifs in 
experimentally determined T-cell epitopes or alternatively using computational 
techniques to predict MHC class 11-btnding peptides and in particular DR-binding 
peptides. 

W098/52976 and WO00/34317 teach computational threading approaches to identifying 
20 polypeptide sequences with the potential to bind a sub-set of human MHC class II DR 
allotypes. In these teachings, predicted T-cell epitopes are removed by the use of 
judicious amino acid substitution within the primary sequence of the therapeutic antibody 
or non-antibody protein of both non-human and human derivation. 

25 Other techniques exploiting soluble complexes of recombinant MHC molecules in 
combination with synthetic peptides and able to bind to T-cell clones from peripheral 
blood samples from human or experimental animal subjects have been used in the art 
[Kern, F. et al (1998) Nature Medicine 4:975-978; Kwok, W.W. et al (2001) TRENDS in 
Immunology 22: 583-588] and may also be exploited in an epitope identification strategy. 



30 



As depicted above and as consequence thereof, it would be desirable to identify and to 
remove or at least to reduce T-cell epitopes from a given in principal therapeutically 
valuable but originally immunogenic peptide, polypeptide or protein. 
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One of these therapeutically valuable molecules is "interleukin-l receptor antagonist (IL- 
IRA)". The protein can be produced by recombinant technologies using a number of 
different host cell types. The amino acid sequence of interleukin-l receptor antagonist 
(IL-IRA) (depicted as one- letter code) is as follows: 

5 RPSGRKSSKMQAFRIWDVNQKTFYLRNNQLVAGYLQGPNVNLEEKIDVVPIE PHALFLGIHGGKMCLSCVK 

SGDETRLQLEAVNITDLSENRKQDKRFAFIRSDSGPTTSFESAACPGWFLCTAMEADOPVSLTNMPDEGVM 
VTKFYFOEDE 

Others have provided IL-IRa molecules [e.g. US 5,075,222], but none have recognized 
the importance of T-cell epitopes to the immunogenic properties of the protein nor have 
10 been conceived to directly influence said properties in a specific and controlled way 
according to the scheme of the present invention. 

However, there is a continued need for interleukin-l receptor antagonist (IL-IRA) 
analogues with enhanced properties. Desired enhancements include alternative schemes 

15 and modalities for the expression and purification of the said therapeutic, but also and 
especially, improvements in the biological properties of the protein. There is a particular 
need for enhancement of the in vivo characteristics when administered to the human 
subject. In this regard, it is highly desired to provide interleukin-l receptor antagonist 
(IL-IRA) with reduced or absent potential to induce an immune response in the human 

20 subject. 

SUMMARY AND DESCRIPTION OF THE INVENTION 

The present invention provides for modified forms of "interleukin-l receptor antagonist 
(IL-IRA)", in which the immune characteristic is modified by means of reduced or 

25 removed numbers of potential T-cell epitopes. The invention discloses sequences 

identified within the IL-IRa primary sequence that are potential T-cell epitopes by virtue 
of MHC class II binding potential. This disclosure specifically pertains a human EL-lRa 
protein being of 152 amino acid residues [Eisenburg, S.P. et al (1991) Proc. Nad Acad. 
Set U.S.A. 88: 5232-5236]. 

30 The invention may be applied to any IL-IRa species of molecule with substantially the 
same primary amino acid sequences as those disclosed herein and would include therefore 
IL-IRa molecules derived by genetic engineering means or other processes and may not 
contain 1 52 amino acid residues. 

IL-IRa proteins such as identified from murine, bovine, canine and other mammalian 
35 sources have in common many of the peptide sequences of the present disclosure and 
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have in common many peptide sequences with substantially the same sequence as those 
of the disclosed listing. Such protein sequences equally therefore fall under the scope of 
the present invention. 

The invention discloses sequences identified within the interleukin- 1 receptor antagonist 
5 (IL- IRA) primary sequence that are potential T-cell epitopes by virtue of MHC class II 

binding potential. This disclosure specifically pertains the human interleukin- 1 receptor 

antagonist (IL-lRA) protein being the 152 amino acid residues. 

The invention discloses also specific positions within the primary sequence of the 

molecule according to the invention which has to be altered by specific amino acid 
10 substitution, addition or deletion without affecting the biological activity in principal. In 

cases in which the loss of immunogenicity can be achieved only by a simultaneous loss of 

biological activity it is possible to restore said activity by further alterations within the 

amino acid sequence of the protein. 

The invention discloses furthermore methods to produce such modified molecules, above 
15 all methods to identify said T-cell epitopes which have to be altered in order to reduce or 
remove immunogenetic sites. 

The protein according to this invention would expect to display an increased circulation 
time within the human subject and would be of particular benefit in chronic or recurring 
disease settings such as is the case for a number of indications for interleukin- 1 receptor 
20 antagonist (IL-lRA). The present invention provides for modified forms of human IL- 
lRA proteins that are expected to display enhanced properties in vivo. These modified IL- 
lRA molecules can be used in pharmaceutical compositions. 



In summary the invention relates to the following issues: 
25 • a modified molecule having the biological activity of interleukin- 1 receptor antagonist 
(IL-lRA) and being substantially non-immunogenic or less immunogenic than any 
non-modified molecule having the same biological activity when used in vivo; 

• an accordingly specified molecule, wherein said loss of immunogenicity is achieved by 
removing one or more T-cell epitopes derived from the originally non-modified 

30 molecule; 

• an accordingly specified molecule, wherein said loss of immunogenicity is achieved by 
reduction in numbers of MHC allotypes able to bind peptides derived from said 
molecule; 

• an accordingly specified molecule, wherein one T-cell epitope is removed; 
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• an accordingly specified molecule, wherein said originally present T-cell epitopes are 
MHC class II ligands or peptide sequences which show the ability to stimulate or bind 
T-cells via presentation on class II; 

• an accordingly specified molecule, wherein said peptide sequences are selected from 
5 the group as depicted in Table I ; 

• an accordingly specified molecule, wherein 1-9 amino acid residues, preferably one 
amino acid residue in any of the originally present T-cell epitopes are altered; 

• an accordingly specified molecule, wherein the alteration of the amino acid residues is 
substitution, addition or deletion of originally present amino acid(s) residue(s) by other 

10 amino acid residue(s) at specific position(s); 

• an accordingly specified molecule, wherein one or more of the amino acid residue 
substitutions are carried out as indicated in Table 2; 

• an accordingly specified molecule, wherein (additionally) one or more of the amino 
acid residue substitutions are carried out as indicated in Table 3 for the reduction in the 

15 number of MHC allotypes able to bind peptides derived from said molecule; 

• an accordingly specified molecule, wherein, if necessary, additionally further alteration 
usually by substitution, addition or deletion of specific amino acid(s) is conducted to 
restore biological activity of said molecule; 

• a DNA sequence or molecule which codes for any of said modified molecule as 
20 specified above and below; 

• a pharmaceutical composition comprising a modified molecule having the biological 
activity of interleukin- 1 receptor antagonist (IL-1RA) as defined above and / or in the 
claims, optionally together with a pharmaceutically acceptable carrier, diluent or 
excipient; 

25 • a method for manufacturing a modified molecule having the biological activity of 
interleukin- 1 receptor antagonist (1L-IRA) as defined in any of the claims of the 
above-cited claims comprising the following steps: (i) determining the amino acid 
sequence of the polypeptide or part thereof; (ii) identifying one or more potential T- 
cell epitopes within the amino acid sequence of the protein by any method including 

30 determination of the binding of the peptides to MHC molecules using in vitro or in 

silico techniques or biological assays; (iii) designing new sequence variants with one 

or more amino acids within the identified potential T-cell epitopes modified in such a 

way to substantially reduce or eliminate the activity of the T-cell epitope as determined 

by the binding of the peptides to MHC molecules using in vitro or in silico techniques 

7 . 
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or biological assays; (iv) constructing such sequence variants by recombinant DN A 
techniques and testing said variants in order to identify one or more variants with 
desirable properties; and (v) optionally repeating steps (ii) - (iv); 
. an accordingly specified method, wherein step (iii) is carried out by substitution, 
> addition or deletion of I - 9 amino acid residues in any of the originally present T-cell 

epitopes; 

. an accordingly specified method, wherein the alteration is made with reference to a 

homologues protein sequence and / or in silica modeling techniques; 
. an accordingly specified method, wherein step (ii) of above is carried out by the 
following steps: (a) selecting a region of the peptide having a known amino acid 
residue sequence; (b) sequentially sampling overlapping amino acid residue segments 
of predetermined uniform size and constituted by at least three amino acid residues 
from the selected region; (c) calculating MHC Class II molecule binding score for each 
said sampled segment by summing assigned values for each hydrophobic amino acid 
residue side chain present in said sampled amino acid residue segment; and (d) 
identifying at least one of said segments suitable for modification, based on the 
calculated MHC Class II molecule binding score for that segment, to change overall 
MHC Class II binding score for the peptide without substantially reducing therapeutic 
utility of the peptide; step (c) is preferably carried out by using a B6hm scoring 
function modified to include 12-6 van der Waal's ligand-protein energy repulsive term 
and ligand conformational energy term by (1) providing a first data base of MHC Class 
II molecule models; (2) providing a second data base of allowed peptide backbones for 
said MHC Class II molecule models; (3) selecting a model from said first database; 
(4) selecting an allowed peptide backbone from said second data base; (5) identifying 
amino acid residue side chains present in each sampled segment; (6) determmmg the 
binding affinity value for all side chains present in each sampled segment; and 
repeating steps (1) through (5) for each said model and each said backbone; 
. a I3mer T-cell epitope peptide having a potential MHC class II binding activity and 
created from immunogenetically non-modified interleukin-l receptor antagonist (IL- 
1RA), selected from the group as depicted in Table I and its use for the manufacture of 
interleukin-l receptor antagonist (IL-IRA) having substantially no or less 
immunogenicity than any non-modified molecule with the same biological activity 

when used in vivo; 
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• a peptide sequence consisting of at least 9 consecutive amino acid residues of a I3mer 
T-cell epitope peptide as specified above and its use for the manufacture of 
interleukin- 1 receptor antagonist (IL- 1 RA) having substantially no or less 

immunogenicity than any non-modified molecule with the same biological activity 
when used in vivo; 

The term "T-cell epitope" means according to the understanding of this invention an 
amino acid sequence which is able to bind MCH EC, able to stimulate T-cells and / or also 
to bind (without necessarily measurably activating) T-cells in complex with MHC II. 
The term "peptide" as used herein and in the appended claims, is a compound that 
includes two or more amino acids. The amino acids are linked together by a peptide bond 
(defined herein below). There are 20 different naturally occurring amino acids involved 
int eh biological production of peptides, and any number of them may be linked in any 
order to form a peptide chain or ring. The naturally occurring amino acids employed in 
the biological production of peptides all have the L-configuration. Synthetic peptides can 
be prepared employing conventional synthetic methods, utilizing L-amino acids, D-amino 
acids, or various combinations of amino acids of the two different configurations. Some 
peptides contain only a few amino acid units. Short peptides, e.g., having less than ten 
amino acid units, are sometimes referred to as "oligopeptides". Other peptides contain a 
large number of amino acid residues, e.g. up to 100 ore more, and are referred to as 
"polypeptides'*. By convention, a "polypeptide" may be considered as any peptide chain 
containing three or more amino acids, whereas a "oligopeptide" is usually considered as a 
particular type of "short" polypeptide. Thus, as used herein, it is understood that any 
reference to a "polypeptide" also includes an oligopeptide. Further, any reference to a 
"peptide" includes polypeptides, oligopeptides, and proteins. Each different arrangement 
of amino acids forms different polypeptides or proteins. The number of polypeptides-and 
hence the number of different proteins-that can be formed is practically unlimited. 
"Alpha carbon (Ca)" is the carbon atom of the carbon-hydrogen (CH) component that is 
in the peptide chain. A "side chain" is a pendant group to Ca that can comprise a simple 
or complex group or moiety, having physical dimensions that can vary significantly 
compared to the dimensions of the peptide. 

The invention may be applied to any interleukin- 1 receptor antagonist (IL- IRA) species 
of molecule with substantially the same primary amino acid sequences as those disclosed 
herein and would include therefore interleukin- 1 receptor antagonist (IL-1RA) molecules 
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derived by genetic engineering means or other processes and may not contain either 1 52 
amino acid residues. 

interleukin-1 receptor antagonist (EL- IRA) proteins such as identified from other 
mammalian sources have in common many of the peptide sequences of the present 
5 disclosure and have in common many peptide sequences with substantially the same 
sequence as those of the disclosed listing. Such protein sequences equally therefore fall 
under the scope of the present invention. 

The invention is conceived to overcome the practical reality that soluble proteins 
10 introduced into autologous organisms can trigger an immune response resulting in 
development of host antibodies that bind to the soluble protein. One example amongst 
others, is interferon alpha 2 to which a proportion of human patients make antibodies 
despite the fact that this protein is produced endogenously [Russo, D. et al (1996) ibid; 
Stein, R. et al (1988) ibid]. It is likely that the same situation pertains to the therapeutic 
15 use of interleukin-l receptor antagonist (EL- IRA) and the present invention seeks to 

address this by providing interleukin-l receptor antagonist (EL- IRA) proteins with altered 
propensity to elicit an immune response on administration to the human host. 

The general method of the present invention leading to the modified interleukin-l 
20 receptor antagonist (IL-IRA) comprises the following steps: 

(a) determining the amino acid sequence of the polypeptide or part thereof; 

(b) identifying one or more potential T-cell epitopes within the amino acid sequence of 
the protein by any method including determination of the binding of the peptides to MHC 
molecules using in vitro or in silico techniques or biological assays; 

25 (c) designing new sequence variants with one or more amino acids within the identified 
potential T-cell epitopes modified in such a way to substantially reduce or eliminate the 
activity of the T-cell epitope as determined by the binding of the peptides to MHC 
molecules using in vitro or in silico techniques or biological assays. Such sequence 
variants are created in such a way to avoid creation of new potential T-cell epitopes by 

30 the sequence variations unless such new potential T-cell epitopes are, in turn, modified in 
such a way to substantially reduce or eliminate the activity of the T-cell epitope; and 
(d) constructing such sequence variants by recombinant DNA techniques and testing said 
variants in order to identify one or more variants with desirable properties according to 
well known recombinant techniques. 

10 
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The identification of potential T-cell epitopes according to step (b) can be carried out 
according to methods describes previously in the prior art. Suitable methods are 
disclosed in WO 98/59244; WO 98/52976; WO 00/343 17 and may preferably be used to 
identify binding propensity of interleukin-1 receptor antagonist (IL-lRA)-derived 
5 peptides to an MHC class II molecule. 



Another very efficacious method for identifying T-cell epitopes by calculation is 
described in the EXAMPLE which is a preferred embodiment according to this invention. 
In practice a number of variant interleukin-l receptor antagonist (IL-IRA) proteins will 
10 be produced and tested for the desired immune and functional characteristic. The variant 
proteins will most preferably be produced by recombinant DNA techniques although 
other procedures including chemical synthesis of interleukin-1 receptor antagonist (IL- 
IRA) fragments may be contemplated. 



15 The results of an analysis according to step (b) of the above scheme and pertaining to the 
human interleukin-1 receptor antagonist (IL-IRA) protein sequence of 152 amino acid 
residues is presented in Table 1 . 

Table 1: Peptide sequences in human interleukin-1 receptor antagonist (IL-IRA) with 
potential human MHC class II binding activity. 

RKSS KMQ A FRIWD, S KMQAFRI W DVNQ , QAFRIW DVNQKT F , 



20 



30 



35 



FRIWDVNQKTFYL, 
WDVNQKTFYLRNN, 
FYLRNNQL VAG YL , 



RIWDVNQKTFYLR 
KT FYLRNNQL VAG 
LRNNQLVAGYLQG 



NQLVAG YLQGPNV, QLVAG YLQG PN VN 

25 AGYLQGPNVNLEE, GYLQGPNVNLEEK 

VNLEEKIDVVPIE, EKIDWPIEPHAL 

DWPI E PH AL FLG , VP I E PH AL FLGI H 

AL FLG I HGGKMCL , L FLG I HGGKMCLS 



GKMCLSCVKSGDE, MCLSCVKSGDETR 



ETRLQLEAVNITD, 
EAVNITDLSENRK, 



TRLQLEAVNITDL 
VNITDLSENRKQD 



ENRKQDKR FAFIR , KR FAFI RS DSGPT 
AFI RS D S G PTT S F , TS FESAAC PG W FL 
PGW FLCT AMEADQ , WFLCTAMEADQPV 



QPVS LTNM P DEG V , 
PDEGVMVTKFYFQ, 



VSLTNMPDEGVMV 

EGVMVTKFYFQEH 
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I W DVNQKT FYLRN, 
TFYLRNNQLVAGY, 
RNNQL VAG YLQG P , 
LVAGYLQGPNVNL/ 
PNVNLEEKI DVVP, 
I DV V PI EPHAL FL , 
HAL FLGI HGG KMC, 
LGIHGGKMCLSCV, 
SCVKSGDETRLQL, 
LQL E AVN I T DLS E , 
TDLS ENRKQDKRF , 
FAFIRSDSGPTTS, 
S FES AACPGW FLC , 
TAMEADQPVSLTN, 
TNMPDEGVMVTKF, 
GVMVTKFYFQEDE 
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Peptides are 13mers, amino acids are identified using single letter code. 



The results of a design and constructs according to step (c) and (d) of the above scheme 
and pertaining to the modified molecule of this invention is presented in Tables 2 and 3. 

Table 2: Substitutions leading to the elimination of potential T-cell epitopes of human 



inter leukin-1 receptor antagonist (IL-IRA) (WT = wild type). 



Residue 
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Residue 
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Table 3: Additional substitutions leading to the removal of a potential T-cell epitope for I 
or more MHC allotypes. 
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Substitution 
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Residue WT 
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Substitution 
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K N 



S T 
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N 



Q R 



Q R 



The invention relates to interleukin- 1 receptor antagonist (IL-1RA) analogues in which 
substitutions of at least one amino acid residue have been made at positions resulting in a 
substantial reduction in activity of or elimination of one or more potential T-cell epitopes 
from the protein. One or more amino acid substitutions at particular points within any of 
the potential MHC class n ligands identified in Table 1 may result in a interleukin-1 
receptor antagonist (IL-1RA) molecule with a reduced immunogenic potential when 
administered as a therapeutic to the human host. Preferably, amino acid substitutions are 
made at appropriate points within the peptide sequence predicted to achieve substantial 
reduction or elimination of the activity of the T-cell epitope. In practice an appropriate 
point will preferably equate to an amino acid residue binding within one of the 
hydrophobic pockets provided within the MHC class II binding groove. 
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It is most preferred to alter binding within the first pocket of the cleft at the so-called P 1 
or PI anchor position of the peptide. The quality of binding interaction between the P 1 
anchor residue of the peptide and the first pocket of the MHC class II binding groove is 
recognized as being a major determinant of overall binding affinity for the whole peptide. 
5 An appropriate substitution at this position of the peptide will be for a residue less readily 
accommodated within the pocket, for example, substitution to a more hydrophilic residue. 
Amino acid residues in the peptide at positions equating to binding within other pocket 
regions within the MHC binding cleft are also considered and fall under the scope of the 
present. 

10 

It is understood that single amino acid substitutions within a given potential T-cell epitope 
are the most preferred route by which the epitope may be eliminated. Combinations of 
substitution within a single epitope may be contemplated and for example can be 
particularly appropriate where individually defined epitopes are in overlap with each 

15 other. Moreover, amino acid substitutions either singly within a given epitope or in 

combination within a single epitope may be made at positions not equating to the "pocket 
residues" with respect to the MHC class II binding groove, but at any point within the 
peptide sequence. Substitutions may be made with reference to an homologues structure 
or structural method produced using in silico techniques known in the art and may be 

20 based on known structural features of the molecule according to this invention. All such 
substitutions fall within the scope of the present invention. 



Amino acid substitutions other than within the peptides identified above may be 
contemplated particularly when made in combination with substitution(s) made within a 

25 listed peptide. For example a change may be contemplated to restore structure or 

biological activity of the variant molecule. Such compensatory changes and changes to 
include deletion or addition of particular amino acid residues from the interieukin- 1 
receptor antagonist (IL-IRA) polypeptide resulting in a variant with desired activity and 
in combination with changes in any of the disclosed peptides fall under the scope of the 

30 present. 



In as far as this invention relates to modified interleukin-l receptor antagonist (IL-IRA), 
compositions containing such modified interleukin-l receptor antagonist (IL-IRA) 
. proteins or fragments of modified interleukin-l receptor antagonist (IL-IRA) proteins and 

15 
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related compositions should be considered within the scope of the invention. In another 
aspect, the present invention relates to nucleic acids encoding modified interleukin- 1 
receptor antagonist (IL-1RA) entities. In a further aspect the present invention relates to 
methods for therapeutic treatment of humans using the modified interleukin- 1 receptor 
5 antagonist (IL- 1 RA) proteins. 



EXAMPLE 

There are a number of factors that play important roles in determining the total structure 
of a protein or polypeptide. First, the peptide bond, i.e., that bond which joins the amino 
10 acids in the chain together, is a covalent bond. This bond is planar in structure, 

essentially a substituted amide. An "amide" is any of a group of organic compounds 
containing the grouping -CONH-. 

The planar peptide bond linking Ca of adjacent amino acids may be represented as 
depicted below: 



15 



25 



\£ :H 

/ 



n: 



Because the OC and the C-N atoms lie in a relatively rigid plane, free rotation does not 
occur about these axes. Hence, a plane schematically depicted by the interrupted line is 
sometimes referred to as an "amide" or "peptide plane" plane wherein lie the oxygen (O), 
carbon (C), nitrogen (N), and hydrogen (H) atoms of the peptide backbone. At opposite 
20 corners of this amide plane are located the Ca atoms. Since there is substantially no 

rotation about the 0=C and C-N atoms in the peptide or amide plane, a polypeptide chain 
thus comprises a series of planar peptide linkages joining the Ca atoms. 
A second factor that plays an important role in defining the total structure or 
conformation of a polypeptide or protein is the angle of rotation of each amide plane 
about the common Ca linkage. The terms "angle of rotation" and "torsion angle" are 
hereinafter regarded as equivalent terms. Assuming that the O, C, N, and H atoms remain 
in the amide plane (which is usually a valid assumption, although there may be some 
slight deviations from planarity of these atoms for some conformations), these angles of 
rotation define the N and R polypeptide's backbone conformation, i.e., the structure as it 
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exists between adjacent residues. These two angles are known as <J> and A set of the 
angles <t>i, where the subscript i represents a particular residue of a polypeptide chain, 
thus effectively defines the polypeptide secondary structure. The conventions used in 
defining the <|>, iy angles, i.e., the reference points at which the amide planes form a zero 
5 degree angle, and the definition of which angle is <J>, and which angle is vj/, for a given 
polypeptide, are defined in the literature. See, e.g„ Ramachandran et al. Adv. Prou Chem. 
23:283-437 (1968), at pages 285-94, which pages are incorporated herein by reference. 
The present method can be applied to any protein, and is based in part upon the discovery 
that in humans the primary Pocket I anchor position of MHC Class II molecule binding 

10 grooves has a well designed specificity for particular amino acid side chains. The 

specificity of this pocket is determined by the identity of the amino acid at position 86 of 
the beta chain of the MHC Class II molecule. This site is located at the bottom of Pocket 
1 and determines the size of the side chain that can be accommodated by this pocket. 
Marshall, K. W., J. Immunol., 152 :4946-4956 (1994). If this residue is a glycine, then all 

15 hydrophobic aliphatic and aromatic amino acids (hydrophobic aliphatics being: valine, 
leucine, isoleucine, methionine and aromatics being: phenylalanine, tyrosine and 
tryptophan) can be accommodated in the pocket, a preference being for the aromatic side 
chains. If this pocket residue is a valine, then the side chain of this amino acid protrudes 
into the pocket and restricts the size of peptide side chains that can be accommodated 

20 such that only hydrophobic aliphatic side chains can be accommodated. Therefore, in an 
amino acid residue sequence, wherever an amino acid with a hydrophobic aliphatic or 
aromatic side chain is found, there is the potential for a MHC Class II restricted T-cell 
epitope to be present. If the side-chain is hydrophobic aliphatic, however, it is 
approximately twice as likely to be associated with a T-cell epitope than an aromatic side 

25 chain (assuming an approximately even distribution of Pocket 1 types throughout the 
global population). 

A computational method embodying the present invention profiles the likelihood of 
peptide regions to contain T-cell epitopes as follows: 

(1) The primary sequence of a peptide segment of predetermined length is scanned, and 
30 all hydrophobic aliphatic and aromatic side chains present are identified. (2)The 

hydrophobic aliphatic side chains are assigned a value greater than that for the aromatic 
side chains; preferably about twice the value assigned to the aromatic side chains, e.g., a 
value of 2 for a hydrophobic aliphatic side chain and a value of 1 for an aromatic side 
chain. (3) The values determined to be oresent are summed for each overlapping amino 
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acid residue segment (window) of predetermined uniform length within the peptide, and 
the total value for a particular segment (window) is assigned to a single amino acid 
residue at an intermediate position of the segment (window), preferably to a residue at 
about the midpoint of the sampled segment (window). This procedure is repeated for 
5 each sampled overlapping amino acid residue segment (window). Thus, each amino acid 
residue of the peptide is assigned a value that relates to the likelihood of a T-cell epitope 
being present in that particular segment (window). (4) The values calculated and assigned 
as described in Step 3, above, can be plotted against the amino acid coordinates of the 
entire amino acid residue sequence being assessed. (5) All portions of the sequence which 
10 have a score of a predetermined value, e.g., a value of I , are deemed likely to contain a T- 
cell epitope and can be modified, if desired. This particular aspect of the present invention 
provides a general method by which the regions of peptides likely to contain T-cell 
epitopes can be described. Modifications to the peptide in these regions have the potential 
to modify the MHC Class II binding characteristics. 
1 5 According to another aspect of the present invention, T-cell epitopes can be predicted 
with greater accuracy by the use of a more sophisticated computational method which 
takes into account the interactions of peptides with models of MHC Class n alleles. The 
computational prediction of T-cell epitopes present within a peptide according to this 
particular aspect contemplates the construction of models of at least 42 MHC Class II 
20 alleles based upon the structures of all known MHC Class II molecules and a method for 
the use of these models in the computational identification of T-cell epitopes, the 
construction of libraries of peptide backbones for each model in order to allow for the 
known variability in relative peptide backbone alpha carbon (Cct) positions, the 
construction of libraries of amino-acid side chain conformations for each backbone dock 
25 with each model for each of the 20 amino-acid alternatives at positions critical for the 
interaction between peptide and MHC Class II molecule, and the use of these libraries of 
backbones and side-chain conformations in conjunction with a scoring function to select 
the optimum backbone and side-chain conformation for a particular peptide docked with a 
particular MHC Class II molecule and the derivation of a binding score from this 
30 interaction. Models of MHC Class II molecules can be derived via homology modeling 
from a number of similar structures found in the Brookhaven Protein Data Bank ("PDB"). 
These may be made by the use of semi-automatic homology modeling software 
(Modeller, Sali A. & Blundell TL.. 1993. J. MolBiol 234:779-815) which incorporates a 
simulated annealing function, in conjunction with the CHARMm force-field for energy 
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minimisation (available from Molecular Simulations Inc., San Diego, Ca.). Alternative 
modeling methods can be utilized as well The present method differs significantly from 
other computational methods which use libraries of experimentally derived binding data 
of each amino-acid alternative at each position in the binding groove for a small set of 
5 MHC Class II molecules (Marshall, K. W., et ai, Biomed. Pept. Proteins Nucleic Acids, 
1(3): 157-162) (1995) or yet other computational methods which use similar experimental 
binding data in order to define the binding characteristics of particular types of binding 
pockets within the groove, again using a relatively small subset of MHC Class II 
molecules, and then 'mixing and matching' pocket types from this pocket library to 
10 artificially create further 'virtual' MHC Class II molecules (Stumiolo T., et al., Nat. 

Biotech, 17(6): 555-561 (1999). Both prior methods suffer the major disadvantage that, 
due to the complexity of the assays and the need to synthesize large numbers of peptide 
variants, only a small number of MHC Class II molecules can be experimentally scanned. 
Therefore the first prior method can only make predictions for a small number of MHC 
15 Class II molecules. The second prior method also makes the assumption that a pocket 
lined with similar amino-acids in one molecule will have the same binding characteristics 
when in the context of a different Class II allele and suffers further disadvantages in that 
only those MHC Class II molecules can be 'virtually' created which contain pockets 
contained within the pocket library. Using the modeling approach described herein, the 
20 structure of any number and type of MHC Class TI molecules can be deduced, therefore 
alleles can be specifically selected to be representative of the global population. In 
addition, the number of MHC Class II molecules scanned can be increased by making 
further models further than having to generate additional data via complex 
experimentation. The use of a backbone library allows for variation in the positions of the 
25 Cot atoms of the various peptides being scanned when docked with particular MHC Class 
II molecules. This is again in contrast to the alternative prior computational methods 
described above which rely on the use of simplified peptide backbones for scanning 
amino-acid binding in particular pockets. These simplified backbones are not likely to be 
representative of backbone conformations found in 'real' peptides leading to inaccuracies 
30 in prediction of peptide binding. The present backbone library is created by superposing 
the backbones of ail peptides bound to MHC Class II molecules found within the Protein 
Data Bank and noting the root mean square (RMS) deviation between the Cot atoms of 
each of the eleven ammo-acids located within the binding groove. While this library can 
be derived from a small number of suitable available mouse and human structures 
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(currently 13), in order to allow for the possibility of even greater variability, the RMS 
figure for each C"-D position is increased by 50%. The average Ca position of each 
ammo-acid is then determined and a sphere drawn around this point whose radius equals 
the RMS deviation at that position plus 50%. This sphere represents all allowed Ca 
5 positions. Working from the Ca with the least RMS deviation (that of the amino-acid in 
Pocket 1 as mentioned above, equivalent to Position 2 of the 1 1 residues in the binding 
groove), the sphere is three-dimensional ly gridded, and each vertex within the grid is then 
used as a possible location for a Ca of that amino-acid. The subsequent amide plane, 
corresponding to the peptide bond to the subsequent amino-acid is grafted onto each of 
10 these Cas and the <\> and \y angles are rotated step-wise at set intervals in order to position 
the subsequent Ca. If the subsequent Ca falls within the 'sphere of allowed positions' 
for this Ca than the orientation of the dipeptide is accepted, whereas if it falls outside the 
sphere then the dipeptide is rejected. This process is then repeated for each of the 
subsequent Ca positions, such that the peptide grows from the Pocket 1 Ca 'seed', until 
15 all nine subsequent Cas have been positioned from all possible permutations of the 
preceding Cas. The process is then repeated once more for the single Ca preceding 
pocket I to create a library of backbone Ca positions located within the binding groove. The 
number of backbones generated is dependent upon several factors: The size of the 
Spheres of allowed positions' ; the fineness of the gridding of the 'primary sphere' at the 
20 Pocket I position; the fineness of the step-wise rotation of the <{> and angles used to 

position subsequent Cos. Using this process, a large library of backbones can be created. 
The larger the backbone library, the more likely it will be that the optimum fit will be 
found for a particular peptide within the binding groove of an MHC Class II molecule. 
Inasmuch as all backbones will not be suitable for docking with all the models of MHC 
25 Class II molecules due to clashes with ainino-acids of the binding domains, for each allele 
a subset of the library is created comprising backbones which can be accommodated by 
that allele. The use of the backbone library, in conjunction with the models of MHC Class 
II molecules creates an exhaustive database consisting of allowed side chain 
conformations for each amino-acid in each position of the binding groove for each MHC 
30 Class II molecule docked with each allowed backbone. This data set is generated using a 
simple steric overlap function where a MHC Class II molecule is docked with a backbone 
and an amino-acid side chain is grafted onto the backbone at the desired position. Each of 
the rotatable bonds of the side chain is rotated step-wise at set intervals and the resultant 
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positions of the atoms dependent upon that bond noted. The interaction of the atom with 
atoms of side-chains of the binding groove is noted and positions are either accepted or 
rejected according to the following criteria: The sum total of the overlap of all atoms so 
far positioned must not exceed a pre-determined value. Thus the stringency of the 
5 conformational search is a function of the interval used in the step-wise rotation of the 
bond and the pre-determined limit for the total overlap. This latter value can be small if it 
is known that a particular pocket is rigid, however the stringency can be relaxed if the 
positions of pocket side-chains are known to be relatively flexible. Thus allowances can 
be made to imitate variations in flexibility within pockets of the binding groove. This 
10 conformational search is then repeated for every amino-acid at every position of each 

backbone when docked with each of the MHC Class II molecules to create the exhaustive 
database of side-chain conformations. A suitable mathematical expression is used to 
estimate the energy of binding between models of MHC Class II molecules in conjunction 
with peptide ligand conformations which have to be empirically derived by scanning the 
15 large database of backbone/side-chain conformations described above. Thus a protein is 
scanned for potential T-cell epitopes by subjecting each possible peptide of length 
varying between 9 and 20 amino-acids (although the length is kept constant for each scan) 
to the following computations: An MHC Class II molecule is selected together with a 
peptide backbone allowed for that molecule and the side-chains corresponding to the 
20 desired peptide sequence are grafted on. Atom identity and interatomic distance data 

relating to a particular side-chain at a particular position on the backbone are collected for 
each allowed conformation of that amino-acid (obtained from the database described 
above). This is repeated for each side-chain along the backbone and peptide scores 
derived using a scoring function. The best score for that backbone is retained and the 
25 process repeated for each allowed backbone for the selected model. The scores from all 
allowed backbones are compared and the highest score is deemed to be the peptide score 
for the desired peptide in that MHC Class II model. This process is then repeated for each 
model with every possible peptide derived from the protein being scanned, and the scores 
for peptides versus models are displayed. 
30 In the context of the present invention, each ligand presented for the binding affinity 
calculation is an amino-acid segment selected from a peptide or protein as discussed 
above. Thus, the ligand is a selected stretch of amino acids about 9 to 20 amino acids in 
length derived from a peptide, polypeptide or protein of known sequence. The terms 
"amino acids" and "residues" are hereinafter regarded as equivalent terms. The ligand, in 

* 
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the form of the consecutive amino acids of the peptide to be examined grafted onto a 
backbone from the backbone library, is positioned in the binding cleft of an MHC Class II 
molecule from the MHC Class II molecule model library via the coordinates of the C- 
□ atoms of the peptide backbone and an allowed conformation for each side-chain is 
5 selected from the database of allowed conformations. The relevant atom identities and 
interatomic distances are also retrieved from this database and used to calculate the 
peptide binding score. Ligands with a high binding affinity for the MHC Class II binding 
pocket are flagged as candidates for site-directed mutagenesis. Amino-acid substitutions 
are made in the flagged ligand (and hence in the protein of interest) which is then retested 
10 using the scoring function in order to determine changes which reduce the binding affinity 
below a predetermined threshold value. These changes can then be incorporated into the 
protein of interest to remove T-cell epitopes. Binding between the peptide ligand and the 
binding groove of MHC Class II molecules involves non-covalent interactions including, 
but not limited to: hydrogen bonds, electrostatic interactions, hydrophobic (lipophilic) 
1 5 interactions and Van der Walls interactions. These are included in the peptide scoring 
function as described in detail below. It should be understood that a hydrogen bond is a 
non-covalent bond which can be formed between polar or charged groups and consists of 
a hydrogen atom shared by two other atoms. The hydrogen of the hydrogen donor has a 
positive charge where the hydrogen acceptor has a partial negative charge. For the 
20 purposes of peptide/protein interactions, hydrogen bond donors may be either nitrogens 
with hydrogen attached or hydrogens attached to oxygen or nitrogen. Hydrogen bond 
acceptor atoms may be oxygens not attached to hydrogen, nitrogens with no hydrogens 
attached and one or two connections, or sulphurs with only one connection. Certain 
atoms, such as oxygens attached to hydrogens or imine nitrogens (e.g. C=NH) may be 
25 both hydrogen acceptors or donors. Hydrogen bond energies range from 3 to 7 Kcal/mol 
and are much stronger than Van der Waal's bonds, but weaker than covalent bonds. 
Hydrogen bonds are also highly directional and are at their strongest when the donor 
atom, hydrogen atom and acceptor atom are co-linear. Electrostatic bonds are formed 
between oppositely charged ion pairs and the strength of the interaction is inversely 
30 proportional to the square of the distance between the atoms according to Coulomb's law. 
The optimal distance between ion pairs is about 2,8 A. In protein/peptide interactions, 
electrostatic bonds may be formed between arginine, histidine or lysine and aspartate or 
glutamate. The strength of the bond will depend upon the pKa of the ionizing group and 
the dielectric constant of the medium although they are approximately similar in strength 
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to hydrogen bonds. Lipophilic interactions are favorable hydrophobic-hydrophobic 
contacts that occur between he protein and peptide ligand. Usually, these will occur 
between hydrophobic amino acid side chains of the peptide buried within the pockets of 
the binding groove such that they are not exposed to solvent. Exposure of the 
5 hydrophobic residues to solvent is highly unfavorable since the surrounding solvent 
molecules are forced to hydrogen bond with each other forming cage-like clathrate 
structures. The resultant decrease in entropy is highly unfavorable. Lipophilic atoms may 
be sulphurs which are neither polar nor hydrogen acceptors and carbon atoms which are 
not polar. Van der Waal's bonds are non-specific forces found between atoms which are 

10 ' 3-4A apart. They are weaker and less specific than hydrogen and electrostatic bonds. 
The distribution of electronic charge around an atom changes with time and, at any 
instant, the charge distribution is not symmetric. This transient asymmetry in electronic 
charge induces a similar asymmetry in neighboring atoms. The resultant attractive forces 
between atoms reaches a maximum at the Van der Waal's contact distance but diminishes 

15 very rapidly at about 1 A to about 2A. Conversely, as atoms become separated by less 
than the contact distance, increasingly strong repulsive forces become dominant as the 
outer electron clouds of the atoms overlap. Although the attractive forces are relatively 
weak compared to electrostatic and hydrogen bonds (about 0.6 Kcal/mol), the repulsive 
forces in particular may be very important in determining whether a peptide ligand may 

20 bind successfully to a protein. In one embodiment, the Bohm scoring function (SCORE 1 
approach) is used to estimate the binding constant. (Bohm, H.J., J. Comput Aided Mol. 
Des. y 8(3):243-256 (1994) which is hereby incorporated in its entirety). In another 
embodiment, the scoring function (SCORE2 approach) is used to estimate the binding 
affinities as an indicator of a ligand containing a T-cell epitope (Bohm, H.J., J. Comput 

25 Aided Mol Des., l2(4):309-323 (1 998) which is hereby incorporated in its entirety). 
However, the Bohm scoring functions as described in the above references are used to 
estimate the binding affinity of a ligand to a protein where it is already known that the 
ligand successfully binds to the protein and the protein/ligand complex has had its 
structure solved, the solved structure being present in the Protein Data Bank ('TDB"). 

30 Therefore, the scoring function has been developed with the benefit of known positive 
binding data. In order to allow for discrimination between positive and negative binders, 
a repulsion term must be added to the equation. In addition, a more satisfactory estimate 
of binding energy is achieved by computing the lipophilic interactions in a pairwise 
manner rather than using the area based energy term of the above Bohm functions. 
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Therefore, in a preferred embodiment, the binding energy is estimated using a modified 
Bolun scoring function. In the modified Bohm scoring function, the binding energy 
between protein and ligand (AG b i n <j) is estimated considering the following parameters: 
The reduction of binding energy due to the overall loss of translational and rotational 

5 entropy of the ligand (AG 0 ); contributions from ideal hydrogen bonds (AGhb) where at 
least one partner is neutral; contributions from unperturbed ionic interactions (AGionic); 
lipophilic interactions between lipophilic ligand atoms and lipophilic acceptor atoms 
(AGiipo); the loss of binding energy due to the freezing of internal degrees of freedom in 
the ligand, i.e., the freedom of rotation about each C-C bond is reduced (AG ro 0; *he 

10 energy of the interaction between the protein and ligand (Evdw). Consideration of these 
terms gives equation 1 : 

(AGticdH AG 0 )+( AGhbXNhtVK AG^^^scM AG^oXNupo)^ AGroc+N^+CE vdw)- 

Where N is the number of qualifying interactions for a specific term and, in one 
embodiment, AG 0 , AG hb , AGionic, AGh po and AG r0 t are constants which are given the 
15 values: 5.4, -4.7, -4.7, -0.17, and 1.4, respectively. 
The term Nhb is calculated according to equation 2 : 

N h b - X* Wsf(AR, Aa) X fCNneighb) X fpcs 

fl(AR, Aa) is a penalty function which accounts for large deviations of hydrogen bonds 
from ideality and is calculated according to equation 3: 
20 f(AR, A-G) = fl(AR) x f2(Aa) 

Where: fl(AR) = 1 if AR<=TOL 

or — 1 - (AR - TOL)/0.4 if AR <=0.4 + TOL 
or =0 if ARXJ.4 + TOL 
And: £2(Aa) -1 ifAa<30° 

25 or - l-( Aa - 30)/50 if Aa <=80° 

or =OifAa>80° 
TOL is the tolerated deviation in hydrogen bond length = 0.25A 
AR is the deviation of the H-O/N hydrogen bond length from the ideal value = 1.9 A 
Aa is the deviation of the hydrogen bond angle Z n/o-h.^ from its idealized value of 
30 180° 

f(N nc ighb) distinguishes between concave and convex parts of a protein surface and 
therefore assigns greater weight to polar interactions found in pockets rather than those 
found at the protein surface. This function is calculated according to equation 4 below: 

fCNncighb) = (N^^ht/Nucigub.0) a where a = 0.5 
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N nc ighb is the number of non-hydrogen protein atoms that are closer than 5 A to any given 

protein atom. 

Noci g hb,o is a constant = 25 

fpc$ is a function which allows for the polar contact surface area per hydrogen bond and 
5 therefore distinguishes between strong and weak hydrogen bonds and its value is 
determined according to the following criteria: 
fpc- B when A^^/Nhs < 10 A 2 
or fpcs— 

1 when Apoa/NHB > 10 A 2 
Apdar is the size of the polar protein- ligand contact surface 
10 N H b is the number of hydrogen bonds 
B is a constant whose value = 1,2 

For the implementation of the modified Bohrn scoring function, the contributions from 
ionic interactions, AGj 0 mc, are computed in a similar fashion to those from hydrogen 
bonds described above since the same geometry dependency is assumed. 
1 5 The term Ni? p <> is calculated according to equation 5 below: 
N lip o = Ztif(ri0 

f(ri L ) is calculated for all lipophilic ligand atoms, 1, and all lipophilic protein atoms, L, 
according to the following criteria: 

f(r» L ) =1 when ri L <= Rlf(r, L ) - Rl)/(R2-Rl) when R2 <r lL > Rl 
20 f(ri L ) =0 when r n , >= R2 

Where: Rl=r, vdw ^r t vdw + 0.5 
and R2= Rl + 3.0 

and ri vdw is the Van der Waal's radius of atom I 
and r L vdw is the Van der Waal's radius of atom L 
25 The term N rot is the number of rotable bonds of the amino acid side chain and is taken to 
be the number of acyclic sp 3 - sp 3 and sp 3 - sp 2 bonds. Rotations of terminal -CH3 or - 
NH 3 are not taken into account. 

The final term, E V dw, is calculated according to equation 6 below: 
Evdw = Ci^r,*" +r 2 v *y 2 /r t2 - (r^ w *t 2 ^V>. where: 
30 ei and £2 are constants dependant upon atom identity 

n vdw +r 2 vdw are the Van der Waal's atomic radii 
r is the distance between a pair of atoms. 

With regard to Equation 6, in one embodiment, the constants 61 and e 2 are given the atom 
values: C: 0.245, N: 0.283, O: 0.316, S: 0.316, respectively (i.e. for atoms of Carbon, 
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Nitrogen, Oxygen and Sulphur, respectively). With regards to equations 5 and 6, the Van 
der Waal's radii are given the atom values C: 1.85, N: 1.75, O: 1.60, S: 2.00A. 
It should be understood that all predetermined values and constants given in the equations 
above are determined within the constraints of current understandings of protein ligand 
5 interactions with particular regard to the type of computation being undertaken herein. 
Therefore, it is possible that, as this scoring function is refined further, these values and 
constants may change hence any suitable numerical value which gives the desired results 
in terms of estimating the binding energy of a protein to a ligand may be used and hence 
fall within the scope of the present invention. As described above, the scoring function is 
10 applied to data extracted from the database of side-chain conformations, atom identities, 
and interatomic distances. For the purposes of the present description, the number of 
MHC Class II molecules included in this database is 42 models plus four solved 
structures. It should be apparent from the above descriptions that the modular nature of 
the construction of the computational method of the present invention means that new 
15 models can simply be added and scanned with the peptide backbone library and side- 
chain conformational search function to create additional data sets which can be 
processed by the peptide scoring function as described above. This allows for the 
repertoire of scanned MHC Class II molecules to easily be increased, or structures and 
associated data to be replaced if data are available to create more accurate models of the 
20 existing alleles. 

The present prediction method can be calibrated against a data set comprising a large 
number of peptides whose affinity for various MHC Class II molecules has previously 
been experimentally determined- By comparison of calculated versus experimental data, 
a cut of value can be determined above which it is known that all experimentally 

25 determined T-cell epitopes are correctly predicted. It should be understood that, although 
the above scoring function is relatively simple compared to some sophisticated 
methodologies that are available, the calculations are performed extremely rapidly. It 
should also be understood that the objective is not to calculate the true binding energy per 
se for each peptide docked in the binding groove of a selected MHC Class II protein. The 

30 underlying objective is to obtain comparative binding energy data as an aid to predicting 
the location of T-cell epitopes based on the primary structure (i.e. amino acid sequence) 
of a selected protein. A relatively high binding energy or a binding energy above a 
selected threshold value would suggest the presence of a T-cell epitope in the ligand. The 
ligand may then be subjected to at least one round of amino-acid substitution and the 
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binding energy recalculated. Due to the rapid nature of the calculations, these 
manipulations of the peptide sequence can be performed interactively within the 

* 

program's user interface on cost-effectively available computer hardware. Major 
investment in computer hardware is thus not required. It would be apparent to one skilled 
5 in the art that other available software could be used for the same purposes. In particular, 
more sophisticated software which is capable of docking ligands into protein binding-sites 
may be used in conjunction with energy minimization. Examples of docking software 
are: DOCK (Kuntz et aL, J. Mol Biol., 161:269-288 (1982)), LUDI (Bohm, H.J., J. 
Comput Aided MoL Des., 8:623-632 (1994)) and FLEXX (Rarey M., etaL. ISMB, 3:300- 

10 308(1995)). Examples of molecular modeling and manipulation software include: 
AMBER (Tripos) and CHARMm (Molecular Simulations Inc.). The use of these 
computational methods would severely limit the throughput of the method of this 
invention due to the lengths of processing time required to make the necessary 
calculations. However, it is feasible that such methods could be used as a Secondary 

15 screen 1 to obtain more accurate calculations of binding energy for peptides which are 
found to be 'positive binders' via the method of the present invention. The limitation of 
processing time for sophisticated molecular mechanic or molecular dynamic calculations 
is one which is defined both by the design of the software which makes these calculations 
and the current technology limitations of computer hardware. It may be anticipated that, 

20 in the future, with the writing of more efficient code and the continuing increases in speed 
of computer processors, it may become feasible to make such calculations within a more 
manageable time-frame. Further information on energy functions applied to 
macromolecules and consideration of the various interactions that take place within a 
folded protein structure can be found in: Brooks, B.R., et al. y J. Comput. Chem., 4:187- 

25 217 (1983) and further information concerning general protein-ligand interactions can be 
found in: Dauber-Osguthorpe et aL, Proteins4(l):3 1-47(1988), which are incorporated 
herein by reference in their entirety. Useful background information can also be found, 
for example, in Fasman, G.O., ed., Prediction of Protein Structure and the Principles of 
Protein Conformation. Plenum Press, New York, ISBN: 0-306 4313-9. 

30 
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Patent Claims 

1 . A modified molecule having the biological activity of inter leukin-l receptor 
antagonist (1L-1RA) and being substantially non-immunogenic or less 
5 immunogenic than any non-modified molecule having the same biological activity 

when used in vivo. 



2. A molecule according to claim U wherein said loss of imrnunogenicity is achieved 
by removing one or more T-celi epitopes derived from the originally non-modified 
10 molecule. 



15 



A molecule according to claim 1 or 2, wherein said loss of imrnunogenicity is 
achieved by reduction in numbers of MHC allotypes able to bind peptides derived 
from said molecule. 



4. A molecule according to claim 2 or 3, wherein one T-cell epitope is removed. 



5. A molecule according to any of the claims 2-4, wherein said originally present T- 
cell epitopes are MHC class II ligands or peptide sequences which show the ability 

20 to stimulate or bind T-cells via presentation on class IL 

6. A molecule according to claim 5, wherein said peptide sequences are selected from 
the group as depicted in Table L 

25 7. A molecule according to any of the claims 2-6, wherein 1 - 9 amino acid residues 

in any of the originally present T-cell epitopes are altered. 

8. A molecule according to claim 7, wherein one amino acid residue is altered. 



30 9. A molecule according to claim 7 or 8, wherein the alteration of the amino acid 

residues is substitution of originally present amino acid(s) residue(s) by other amino 
acid residue(s) at specific posit ion(s). 
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10. A molecule according to claim 9, wherein one or more of the amino acid residue 
substitutions are carried out as indicated in Table 2. 



11. A molecule according to claim 10, wherein additionally one or more of the amino 
5 acid residue substitutions are carried out as indicated in Table 3 for the reduction in 

the number of MHC allotypes able to bind peptides derived from said molecule. 



12. A molecule according to claim 9, wherein one or more amino acid substitutions are 
carried as indicated in Table 3. 



10 



13. A molecule according to claim 7 or 8, wherein the alteration of the amino acid 
residues is deletion of originally present amino acid(s) residue(s) at specific 
position(s). 

15 14. A molecule according to claim 7 or 8, wherein the alteration of the amino acid 
residues is addition of amino acid(s) at specific position(s) to those originally 
present. 

15. A molecule according to any of the claims 7 to 14, wherein additionally further 
20 alteration is conducted to restore biological activity of said molecule. 

s 

16. A molecule according to claim 15, wherein the additional further alteration is 
substitution, addition or deletion of specific amino acid(s). 

25 17. A modified molecule according to any of the claims 7-16, wherein the amino acid 

alteration is made with reference to an homologous protein sequence. 

18. A modified molecule according to any of the claims 7-16, wherein the amino acid 
alteration is made with reference to in silico modeling techniques. 



30 



19. A DNA sequence coding for a modified interleukin-l receptor antagonist (IL-1RA) 
of any of the claims 1-18. 
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20. A pharmaceutical composition comprising a modified molecule having the 
biological activity of interleukin-1 receptor antagonist (IL-IRA) as defined in any 
of the above-cited claims, optionally together with a pharmaceutical^ acceptable 
carrier, diluent or excipienL 

21. A method for manufacturing a modified molecule having the biological activity of 
interleukin-1 receptor antagonist (IL-IRA) as defined in any of the claims of the 
above-cited claims comprising the following steps: 

(i) determining the amino acid sequence of the polypeptide or part thereof 

(ii) identifying one or more potential T-cell epitopes within the amino acid sequence 
of the protein by any method including determination of the binding of the peptides 
to MHC molecules using in vitro or in silico techniques or biological assays; 

(iii) designing new sequence variants with one or more amino acids within the 
identified potential T-cell epitopes modified in such a way to substantially reduce or 
eliminate the activity of the T-cetl epitope as determined by the binding of the 
peptides to MHC molecules using in vitro or in silico techniques or biological 
assays, or by binding of peptide-MHC complexes to T-cells; 

(iv) constructing such sequence variants by recombinant DNA techniques and 
testing said variants in order to identify one or more variants with desirable 
properties; and 

(v) optionally repeating steps (ii) - (iv). 

22. A method of claim 21 , wherein step (iii) is carried out by substitution, addition or 
deletion of 1 - 9 amino acid residues in any of the originally present T-cell epitopes. 

23. A method of claim 22, wherein the alteration is made with reference to a 
homologues protein sequence and / or in silico modeling techniques. 

24. A method of any of the claims 21-23, wherein step (ii) is carried out by the 
following steps: (a) selecting a region of the peptide having a known amino acid 
residue sequence; (b) sequentially sampling overlapping amino acid residue 
segments of predetermined uniform size and constituted by at least three amino acid 
residues from the selected region; (c) calculating MHC Class II molecule binding 
score for each said sampled segment by summing assigned values for each 
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hydrophobic amino acid residue side chain present in said sampled amino acid 
residue segment; and (d) identifying at least one of said segments suitable for 
modification, based on the calculated MHC Class II molecule binding score for that 
segment, to change overall MHC Class II binding score for the peptide without 
5 substantially reducing therapeutic utility of the peptide. 

25. A method of claim 24, wherein step (c) is carried out by using a Bohm scoring 

function modified to include 12-6 van der Waal's ligand-protein energy repulsive 
term and ligand conformational energy term by (1) providing a first data base of 

10 MHC Class II molecule models; (2) providing a second data base of allowed 

peptide backbones for said MHC Class II molecule models; (3) selecting a model 
from said first data base; (4) selecting an allowed peptide backbone from said 
second data base; (5) identifying amino acid residue side chains present in each 
sampled segment; (6) determining the binding affinity value for all side chains 

15 present in each sampled segment; and repeating steps (1) through (5) for each said 

model and each said backbone. 



26. A 13mer T-cell epitope peptide having a potential MHC class II binding activity 
and created from immunogenically non-modified interleukin-l receptor antagonist 
20 (IL-1RA), selected from the group as depicted in Table 1 . 



27. A peptide sequence consisting of at least 9 consecutive amino acid residues of a 
I3mer T-cell epitope peptide according to claim 26. 

25 28. Use of a I3mer T-cell epitope peptide according to claim 26 for the manufacture of 

interleukin-l receptor antagonist (IL-1RA) having substantially no or less 
immunogenicity than any non-modified molecule with the same biological activity 
when used in vivo. 



30 29. Use of a peptide sequence according to claim 27 for the manufacture of interleukin- 

1 receptor antagonist (IL-1RA) having substantially no or less immunogenicity than 
any non-modified molecule with the same biological activity when used in vivo. 
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